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A mechanism is proposed based on the Kubo formula to generate a spin-polarized magneto-optical 
current of Dirac electrons in solids which have strong spin-orbit interactions such as bismuth. The ac 
current response functions are calculated in the isotropic Wollf model under an external magnetic field, 
and the selection rules for Dirac electrons are obtained. By using the circularly polarized light and tuning 
its frequency, one can excite electrons concentrated in the spin-polarized lowest Landau level when the 
chemical potential locates in the band gap, so that spin-polarization in the magneto-optical current can be 
achieved. 
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The motion of electrons can be described by a Dirac type 
equation when there are two bands coupled by a strong spin- 
orbit interaction as shown by Wolff in the study of bismuth 
and its alloys.'"^' The validity of this approach has been con- 
firmed experimentally. In the Wolff Hamiltonian (three- 
dimensional 4x4 matrix and massive), the spin-orbit interac- 
tion is included in the off-diagonal elements as in the original 
Dirac Hamiltonian in relativistic quantum mechanics. These 
spin-dependent off-diagonal elements lead to non-trivial inter- 
band effects such as large diamagnetism^' and unconven- 
tional inter-band Hall effects.^' Although the graphene-type 
Hamiltonian (e.g., in graphene' or Q'-ET2l3'^*) also gives 
Dirac electrons, they are intrinsically different from the Wolff 
Hamiltonian with respect to the spin properties. The graphene 
Hamiltonian (two-dimensional 2x2 matrix and massless) is 
not related to the real spin, but to the pseudo spin, which ex- 
presses the degrees of freedom of the sublattice. It has nothing 
to do with spin-orbit interactions or spin properties. Wolff dis- 
cussed in his Hamiltonian that the spin-orbit interaction gives 
rise to photo-induced spin transitions in magneto-optical mea- 
surements. This spin transition itself was an important new 
feature at that time, but no particular spin phenomena have 
been proposed so far using the spin-orbit interaction or the 
spin transition. 

In the present Letter, we propose a mechanism to generate 
a spin-polarized electric current by use of the lowest Landau 
level and the circularly polarized light. Before showing the 
details of calculations, we illustrate the phenomena in terms 
of the Landau levels (Fig. 1). In the presence of strong spin- 
orbit interaction as in bismuth, namely, in Wolff Hamiltonian, 
the energy dispersion when there is no magnetic field (B = 0), 
and the Landau levels in the case of B > are shown in Fig. 1 
(a). The Landau levels are given by Ej = + -y/A^ + 2Aa>*j, 
where 2A is the band gap and w* is the renormalized cy- 
clotron frequency (origin of energy is taken at the center of 
the band gap).''^' The interval of Landau levels decreases as 
j increases. The total angular momentum j is related to the 
orbital- and spin-quantum numbers as j - n -H 1/2 + cr/2, 
where a - -i-l and - 1 correspond to the up and down spin of 



(a) bismuth type 



(b) graphene type 
(for comparison) 






Fig. 1. (color online) (a) Illustrations of Landau levels for bismuth. The 
numbers beside the levels denoting the total angular momentum j, the + 
and - correspond to cr, and "c" and "v" are the abbreviations of conduction 
and valence band, respectively, (b) (for comparison) The Landau levels for 
graphene, which are given by E„ = sgn(«)y V^eSNTt', where y is the 
velocity.'' There, double degeneracy is due to the degrees of freedom of 
sublattice pseudo spin, not the real spin. The possible transitions from the 
selection rule are displayed by an'ows for (c) insulating (|/i| < A), and (d) 
metallic (ji/A = 2.5) regions. 



electrons. For j > 1, the Landau levels with n and cr = -1, 
and with n - 1 and cr = h-1 are degenerate. On the other hand, 
the lowest Landau level with 7 = permits only the state with 
n - and cr - -1. This lowest Landau level with a spin 
of particular direction is a striking characteristic of the Wolff 
Hamiltonian. Such a state cannot be obtained in the case of 
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gapless graphene Hamiltonian as in Fig. 1 (b). 

We discuss optical conductivities in the Landau levels 
shown in Fig. 1 (a). Since the dipole matrix elements are rele- 
vant to the optical conductivity, it turns out that the transitions 
of AJ = + 1 are allowed in the optical conductivity, o-xx{(^), as 
is depicted in the left panels of Figs. 1 (c) and (d). In the in- 
sulating state (Fig. 1 (c)), we see that there are two processes 
for the smallest excitation energy involving 7 = states. In 
this case we do not have net spin polarization. On the other 
hand, if we use a circularly polarized light, i.e., in the case 
of cr+-{(jS), the allowed excitations are restricted to that with 
Ay = -I by the selection rule. Therefore, for the smallest ex- 
citation energy, only the transition shown in the right panel of 
Fig. 1 (c) is allowed, which leads to the spin polarization. 

For our spin-polarized electric current, the spin-orbit inter- 
action is essential (actually, the spin transition always accom- 
panies the orbital transition), and the Landau levels are char- 
acterized by j not by n. Our mechanism uses (I) the inter- 
band transition, (2) the lowest Landau level, and (3) the cir- 
cularly polarized light, whereas Wolff mainly investigated the 
intra-band transition and did not consider the specific prop- 
erty of the lowest Landau level. In graphene, in contrast, it 
is not possible to obtain the spin polarization by the present 
mechanism, since its Hamiltonian does not have spin-orbit in- 
teractions or spin transitions. Note also that the present spin- 
polarized current is (i) a current going through the bulk and 
not along the surface, (ii) accompanied by an electric current 
(i.e., not a pure spin current), and (iii) realized under the mag- 
netic field. These properties are different from the spin Hall 
effects,'^' "'^ the quantum spin Hall effects or the surface state 
of the topological insulators.'^' 

In the following, we show explicit calculations based on 
the "isotropic" Wolff Hamiltonian. As an effective model of 
bismuth, Cohen and Blount employed a simple one electron 
Hamiltonian with the spin-orbit interaction:^' 

.2 v2y 



2m Sm-^c-^ Arrfic^^ ^ ^ 



(1) 



The Pauli matrix cr corresponds to the real spin of electrons, 
V is the crystal potential, and the last term expresses the spin- 
orbit interaction. Wolff found that eq. (1) can be written in 
an essentially identical form to the 4x4 Dirac Hamiltonian 
on the basis of the k • p theory.'' The velocity of bismuth is 
anisotropic by nature. However in this Letter, we discuss the 
case with isotropic velocity as 



A \yk ■ cr 
-\yk ■ cr -A 



(2) 



which we call the isotropic Wolff Hamiltonian. The wave vec- 
tor k is measured from k^, where the bands have their ex- 
trema. (Such a situation is realized at the L-point in bismuth.) 
The spin-orbit interaction is included in the off-diagonal el- 
ements and its magnitude is proportional to y. The paraboUc 
part 12m is negligibly small, when the spin-orbit interaction 
is much larger than the band gap (this is true for bismuth). 

Under an external magnetic field, the eigen energy can be 
easily obtained by the squared wave equation: 



A^+2A(=^ 







(3) 



Here the 2 x 2 Hamiltonian, M", 



2nC 



(4) 



is the same as that for a free particle in a magnetic field (tt - 
p + eA/c) with the cyclotron mass m* = A/y^ and the g-f actor 
g* - Imlml. (jiB = e/2mc is the Bohr magneton, and e > 0.) 
The eigen energy of * is exactly given in the well-known 
manner, and then we have 

±E„,^ik,) = ± ^A2 + 2A [{n + (1 + cr)/2K + kl/2ml], (5) 

where the magnetic field is along the z-direction, a>l = 
eB/mlc, and cr = +1 is the eigenvalues of cr^. The plus and 
minus sign in eq. (5) correspond to the conduction and va- 
lence band, respectively. The energy levels given by eq. (5) 
with A:^ = are shown in Fig. 1 (a). 

It should be emphasized here that the conduction and va- 
lence band have different B-dependences in eq. (5). In the 
conduction band, the energy for cr = -hI (up spin) increases, 
while that for cr = -1 (down spin) decreases (see Fig. 1 (a)). 
In the valence band, on the other hand, the energy for cr = -i-l 
decreases, while that for cr = - 1 increases. In other words, the 
sign of the g-factor is opposite between the conduction and 
valence bands. Furthermore, the spin and cyclotron masses 
are the same, so that the spin splitting is equal to the orbital 
sphtting and becomes very large since m* is small. Then, the 
lowest Landau level is clearly separated from the opposite 
spin levels even by a very weak field. 

The magneto-optical conductivity tensor is given on the ba- 
sis of the Kubo formula in the form'*' 



(Tuvioj) 
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.(o) + i6) - %v(0 + i<J)] , 



(6) 



%yiioJd = -e^T Yj{i\v^\j){j\vy\i)msn)^jii£n - iojd, (7) 

where v = dJif/dk is the velocity operator, ^^(ie„) = 
{i\ [ie„ - ]"' 10 is the Green function, and \i} is the state 
in which the energy is given by Ei{kJ. The effect of damping 
is expressed by F as e„ = e„ + sgn(e„)r. For the inter-band 
(valence — > conduction) transition, we have 



F(aj,-En,o;En+i,o-)iA 



n+lo",HO"- 



X nuolin + 1) [(£„+i,cr - £„,cr + 2A) -i- (r(£„+i,^ + £„,^)]^ 

+ F(a>, £„_i,,)(A-i^_„^)^ 

X mcj'n [(E„-\ + 2A) - (r{En-l,<r + En,a)] 



+ 2F(w, £„,-^)(A^! 



(8) 



where A^l = eBI2nc is the degrees of the Landau-level degen- 
eracy, A^;,,,^ = [En^En'AEna + A)(£„v' + A)]'"^ , and £„,^ 
implicitly includes the dependences as En,^{k^. The energy 
—En,a- corresponds to the initial state in the valence band, and 
En±i,cr and £„ -o- to the excited state in the conduction band. 
We performed the analytic continuation of the Green function 
part 5f(ie„)^^(ie„ - iw^), and obtained the analytical form at 
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(a) Insulating (|x / A = 0.0) 
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Fig. 2. (color online) Frequency dependence of Re o",., (left) and Re o"+ 
(right) for (a) p/A = 2.5 and (b) /i/A = 0.0. 
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For the intra-band (conduction band) transition, we have the 
resuh just by replacing -E„^cr — > +E„,a- in eq. (8). In eq. (8), 
each term expresses the transition -E„cr — > E„-t[^cr (orbital 
transition) and -E„,j- — > £„ _o- (spin transition), which are all 
Aj - ±1 transitions; this is the selection rule for cr_xx- (The 
quantum numbers n and cr are no more good quantum num- 
bers, but only j is the good quantum number) 

The results of crvv(w) are shown in the left of Fig. 2 
for (a) insulating (yu - 0.0) and (b) metallic (/^/A = 2.5) 
cases. All calculations are carried out for the whole processes 

(O':': + + ^cv + (j,vc-) ^jjj^ ^ ^ j q ^jjJ p/^ ^ q q2 

Each peak structure can be understood from the selection rule 
above. For the insulating region, fi/A - 0.0, all peaks appear 
in the high energy (w > 2A) region, which are the contribu- 
tion from the inter-band transitions of (j - Oy — > Ic) and 
ij = Iv ^ Oc) at ll>u (j = Iv -> 2c) and (/ = 2v -> Ic) at W2, 
etc. (Fig. 1 (c)). (The subscripts denote the conduction (c) and 
valence (v) band levels, and cjj = Ej-\{k, = 0) -i- Ej{k, - 0) 
as is indicated by arrows in Fig. 2.) For the metallic region, 
yu/A = 2.5, the high energy peaks are due to the inter-band 
transitions of (j = 2v — » 3c) at ciJt,, (j = 3v — > 4c) and 
ij = 4v — > 3c) at 0)4, etc. (Fig. 1 (d), fi locates between 
j - 2 and 3 in the conduction band). The low energy {a> < A) 
sharp peak is the contribution from the intra-band transition 
of (j = 2c — » 3c). In both regions, the excited states consist of 
both spins at every ajj, so that the current is not spin-polarized. 
Therefore, in crvx(aj), remarkable spin polarization cannot be 
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Fig. 3. (color online) Frequency dependence of Re as^^ (left) and Re crs+ 
(right) for /i/A = 0.0. 



obtained. However, the situation completely changes if we see 
the response to the circularly polarized light. 

The response to the circularly polarized light can be cal- 
culated by the same manner as in (t^x in terms of Vj- = 
(vx + ivy)/ V2. The inter-band contribution to <t>+- is then ob- 
tained as 

^ F(aj,-E„,^,E„^,M''-i.,nJ 



n,k',o 



X mu>ln [{E„-\^a - En,cr + 2A) - Cr(En-\.cr + En,cr)f 



+ IFioj, £„,_i)(A^^_i.„^i)\?(£„,+, + En,-x) 



(10) 



In this case, we have only the transitions of Ay - -1. The re- 
sults of cr+_(w) are shown in the right of Fig. 2. At first sight, 
there is no essential difference between crxx{<jS) and cr+^iu)), 
but their spin structures are definitely different. If we see 
the lowest excited state at wi in the insulating region, only 
ij = Iv — > Oc) transition is possible for cr^_, while both 
ij = Ov — > Ic) and ij = Iv — > Oc) are possible for cTxx 
(Fig. 1 (c)). All electrons excited by circularly polarized light 
of 10 - (jL)\ are concentrated at the lowest Landau level with 
a particular spin direction, i.e., the spin is polarized. Such a 
transition can be realized due to the spin transition in addition 
to the orbital transition. If we have orbital transition only, both 
spins are excited into the j - Oc level by wi, since the initial 
state j - Ov consists of both spins, and their spins are con- 
served only by the orbital transition. In the present case, how- 
ever, the spin-transition is also possible due to the spin-orbit 
coupling, so that the orbital (1^ O^T ) and spin (0^ — > O^T) 
transition lead the excited state consisting of down spin only. 
Therefore, we can obtain the spin-polarized current by the cir- 
cularly polarized light. 

In order to verify the spin polarization of cr+-ico), let us cal- 
culate the current of the magnetic moment, which originates 
from the spin of excited electrons. As discussed before, the 
sign of g-factor in the conduction and valence band are oppo- 
site, so that we define the magnetic moment /ie as 



A'e = -g m 





-(T 



(11) 



In addition, we define the velocity of the magnetic moment 
along the z-direction by Vs/ - iig,Vilg*HQ [i = x,y), which is 
an Hermitian operator. (Here, v, being the velocity operator 
introduced in eq. (7)) Then, the spin-polarization of the cur- 
rent can be estimated with this velocity of the magnetic mo- 
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ment as (D^^v = -e^T En,;,/'>'s//l j><i1vv|/)^(ie„)f#/ie„ - \(jJa)- 
It should be noted that the spin current is not well defined 
since the spin density is not a conserved quantity in the pres- 
ence of spin-orbit coupling. Here we use Os^y just to see the 
tendency of the spin-polarization. The responses 3)s«, ,SH — 
calculated in the same manner as ^xx,+-, and their results are 
shown in Fig. 3 for the insulating region. 

The spin polarization is evaluated by Ics^yl/cr^y. Therefore, 
Figs. 2 and 3 signify that cr+_(w) in the insulating region is 
highly spin-polarized, while cr„(w) is not. Especially, at the 
lowest excitation energy, |crs+-(wi)|/cr+_(a>i) = 1, namely, we 
have 100% spin polarization as is expected. (More precisely, 
the spin-polarization becomes slightly less than 100% for fi- 
nite F, e.g., |crs+-|/cr+_ = 0.97 at wi for F/A = 0.02.) 

One might think that the spin relaxation is very fast in cases 
with a strong spin-orbit coupling, thereby destroying the spin- 
polarization, but this is not the case here. The Landau lev- 
els under consideration are the pure eigenstates of the Wolff 
Hamiltonian, which already includes the effect of strong spin- 
orbit coupling, so that different Landau levels are not mixed 
with each other, i.e., the quantum number j conserves. More- 
over, the present spin-polarized current uses the lowest Lan- 
dau level, whose spin is unique and is not mixed, so that the 
spin-polarization is kept. 

Now we discuss the implications of the present results on 
the experimental measurements. Our spin-polarized current 
flows through the bulk not the surface current as in the topo- 
logical insulators. '^^ The magnitude of the current is large, 
and it should be much easier to detect. Actually, the magneto- 
optical measurements on bismuth already exhibit clear peak 
structures.^* Therefore, the present proposal of spin-polarized 
current will be confirmed experimentally, e.g., by the ordi- 
nary magneto-optical measurements with the use of the Ken- 
rotation microscopy."** There, we will see the bulk spin polar- 
ization through the photo-induced change in magnetization. 

The spin-orbit interaction of bismuth (~ 1.5 eV''*) is the 
largest among the non-radioactive elements, and much larger 
than the band gap 2A ^ 15 meV. This is the ideal situation 
of the Wolff Hamiltonian. The cyclotron mass is extremely 
small m* ^ 0.002 (B \\ bisectrix axes''*), and the g-factor ex- 
ceeds 1000, which enables the huge spin splittings. The sit- 
uation in Fig. 2 (a>* = A) is achieved only by B = 0.12 
T. The insulating Dirac electrons appear in Bii_jSbj alloys 
(0.07 < X < 0.22).2"-22) Also, CaaPbO and its family would 
be good candidates for the insulating Dirac electron.-^^* 

We discussed the transport phenomena with the isotropic 
Wolff Hamiltonian in this Letter. However, the velocity of bis- 
muth is highly anisotropic, so that it would be important to 
take into account the anisotropy when we compare with ex- 
perimental results.^' ^' Generally speaking, if we consider the 
anisotropy, other matrix elements can appear, which will give 
additional orbital- and spin-transition terms. The clarification 
of these consequences needs further studies. 

The Wolff Hamiltonian is essentially different from the 
graphene Hamiltonian with respect to the spin properties. In 
the Wolff Hamiltonian, the signs of g-factor for the real spin 
of electrons are opposite between the conduction and valence 
band, while they are the same in the graphene Hamiltonian. 
Furthermore, the spin splitting is extremely large and the low- 
est Landau level includes only a single spin with a particular 



direction in the former, while the spin splitting is tiny and the 
lowest Landau level is spin degenerate in the latter Therefore, 
the spin-polarized current discussed in this Letter is only pos- 
sible in the Wolff Hamiltonian, while it is difficult to realize in 
the graphene Hamiltonian. Nevertheless, the spin properties 
of the Wolff Hamiltonian are closely related to the pseudo- 
spin properties of the graphene Hamiltonian,^'*'^^* so that an 
interesting pseudo-spin-polarized current may be reahzed by 
the present mechanism in the graphene Hamiltonian. 

We have predicted the spin-polarized electric current in the 
isotropic Wolff model by calculating o-xx{<j^) and cr+_(w) on 
the basis of the Kubo formula. It has been shown that the spin- 
polarized current is possible by the circularly polarized light 
in the insulating regions. The keys to achieving the spin polar- 
ization are (i) the well-isolated lowest Landau level originated 
from the large spin splitting; (ii) the non-trivial spin transi- 
tion in addition to the usual orbital transition; and (iii) the op- 
posite sign of the g-factor or the magnetic moment between 
the conduction and valence band. These points are originated 
from the strong spin-orbit interaction, and are the prominent 
characteristics of the Dirac electrons with spin-orbit interac- 
tions. It would be exceedingly interesting to test these ideas 
by magneto-optical measurements with the use of Kerr rota- 
tion microscopy. 

We thank J. Inoue, K. Miyake and H. Kohno for helpful 
comments. This work is financially supported by MRL Sys- 
tem, and by "HISHO" The Top Thirty Young Researchers of 
Osaka University Project, and also by Young Scientists (B) 
(No. 23740269) from Japan Society for the Promotion of Sci- 
ence. 



1) p. A. Wolff: J. Phys. Chem. Sohds 25 (1964) 1057. 

2) M. H. Cohen and E. I. Blount: Phil. Mag. 5 (1960) 115. 

3) M. S. Dresselhaus: J. Phys. Chem. Solids 32 (1971) 3. 

4) V. S. Edel'man: Adv. Phys. 25 (1976) 555. 

5) M. P. Vecchi, J. R. Pereira, and M. S. Dresselhaus: Phys. Rev. B 14 
(1976) 298. 

6) Z. Zhu, B. Fauque, Y. Fuseya, and K. Behnia: arXiv:l 107.25 17vl . 

7) H. Fukuyama and R. Kubo: J. Phys. Soc. Jpn. 28 (1970) 570. 

8) Y. Fuseya, M. Ogata, and H. Fukuyama: Phys. Rev. Lett. 102 (2009) 
066601. 

9) Y. Zheng and T. Ando: Phys. Rev. B 65 (2002) 245420. 

10) V. P Gusynin and S. G. Sharapov: Phys. Rev. B 73 (2006) 245411. 

11) H. Fukuyama: J. Phys. Soc. Jpn. 76 (2007) 043711. 

12) A. Kobayashi, Y. Suzumura, and H. Fukuyama: J. Phys. Soc. Jpn. 77 
(2008) 064718. 

13) S. Murakami, N. Nagaosa, and S. C. Zhang: Science 301 (2003) 1348. 

14) J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwrrth, and A. H. 
MacDonald: Phys. Rev. Lett. 92 (2004) 126603. 

15) C. L. Kane and E. J. Mele: Phys. Rev. Lett. 95 (2005) 146802. 

16) H. Fukuyama: Prog. Theor Phys. 42 (1969) 1284. 

17) H. Zhang, C.-X. Liu, Z.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang: Nature 
Phys. 5 (2009) 438. 

18) Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. Awschalom: Science 
306 (2007) 1910. 

19) Y. Liu and R. E. Allen: Phys. Rev. B 52 (1995) 1566. 

20) L. Wehrii: Phys. Kondens. Materie 8 (1968) 87. 

21) N. B. Brandt, M. V. Semenov, and L. A. Falkovsky: J. Low. Temp. Phys. 
27 (1977) 75. 

22) B. Lenoir, M. Cassart, J. -P. Michenaud, H. Scherrer, and S. Scherrer: J. 
Phys. Chem. Solids 57 (1996) 89. 

23) T. Kariyado and M. Ogata: J. Phys. Soc. Jpn. 80 (201 1) 083704. 

24) M. Koshino and T. Ando: Phys. Rev. B 81 (2010) 195431. 

25) M. Koshino and T. Ando: Solid State Comm. 151 (2011) 1054. 



